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Abstract 
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1 Introduction 

As is observed [1] , there has been in the last decade active interest to application of fractional 
derivatives in various fields of physics, where they play essential and sometimes leading role in 
understanding the complex classical and quantum systems. Surprisingly, the achievements of the 
fractional-differential calculus in the theory of gravitation and cosmology are very insufficient and 
probably no more then thirty or so publications are devoted to this subject (see, for example, 
[2]- [7]). At the same time, some essential problems in modern theoretical cosmology [8, 9] are 
known that lead to serious modification of the possible sources of the accelerated expansion of 
the Universe (or the gravitational theory itself). Under these circumstances, the interest to the 
■ fractional-differential modifications of cosmological models that appeared recently is fully justified. 

In this paper, we provide a review of the main publications devoted to the fractional differential 
CO ■ approach in cosmology (known to the author). We also offer a novel approach to obtaining the 

modified Friedmann equations consistent with main principles and allowing, in our opinion, to 
avoid the conflict between the integer dimension of (pseudo) Riemann space - time of GR and 
fractional order of derivatives in the modified equations. As is mentioned in [5], there are two 
different methods of approaching what fractional derivative cosmology could be. The simplest 
is the Last Step Modification (LSM) method, in which the Einstein's field equations for a given 
configuration are replaced with analogous fractional field equations. In other words, dt —> Df after 
the field equations for a specific geometry have been derived. The fundamentalist methodology 
is the First Step Modification (FSM), in which one starts by constructing fractional derivative 
geometry. The intensively developed approach to modification of the main cosmological equations 
and non-conservative systems of Lagrangian dynamics on the basis of a variational principle for 
the action of a fractional order (Fractional Action-Like Variational Approach FALVA) developed 
in [2]- [4] represents one of the possible version of intermediate modification (Intermediate Step 
Approach, ISA) mentioned in [5]. 

Our work is based in its main part on ISA and FALVA, because the equations of standard 
cosmology can be obtained from a variational principle for the Einstein-Hilbert action, in which 
the variation is made both over the scale factor and a laps function N, or, as is shown in [10], 
the last can be replaced by a condition of a time scale invariancy, and in any case it is easy 
to generalize masters equations to the case of fractional derivatives. We obtain the fractional- 
differential analogue of the Friedmann equations including a scalar field as a source of gravitation 
in frameworks of ISA and on the basis of FALVA. Until now in cosmology with fractional derivatives, 
there exist more questions than answers to such of them as: How should such models be obtained? 
How should the basic equations of these models be written down? What is the meaning of these 
solutions for the modified models in aspect of the modern problems of cosmology? In our work, 
we try to answer some of these questions. 

1 E-mail: vkshch@yahoo.com 
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2 Fractional integrals and derivatives 



There are no less than two dozens of definitions of the fractional derivative [1]. In the physical 
applications of the fractional differential calculus, more often one deals with Riemann-Liouville 
derivative (RLD), Caputo derivative (CD), and some others. 

Such derivatives are defined by means of analytical continuation of the Cauchy formula for the 
multiple integral of integer order as a single integral with a power-law core into the field of real 
order ^ > 0: 



X 



a) 



The Riemann-Liouville derivative of fractional order a > of function f{x) is defined as the integer 
order derivative of the fractional-order integral (1): 



J JP_ f fit) 

T(n - a) dx n J (x - t) c 



D°f(x) ee Dl c I^f{x) = ^—3— / ,„ ^' dt (2) 



where D™ = d n /dx n , n = [a] + 1. This definition corresponds to the so-called left derivative, 
frequently denoted as c D"f(x). For the limit a = 1, this definition gives df(x)/dx. For example, 
the left RLD of x k for a < 1, c = equals: 

n a k _ r ( fc + 1 ) T k- a l o\ 

U * X ~T(k + l-a) X ■ {6) 

For a = 1, one has the usual result: D\x k — kx k ~ 1 . The interesting feature of RLD is that RLD 
of non-zero constant Co does not equal zero, but for a < 1 it equals -D"Co = Cox~ a /T(l — a). 
The right RLD is defined similarly to (2) on the interval [c, d]: 

d 

x 

The Caputo derivative is the integral transform of a regular derivative, and it is defined by 
moving the integer-order derivative in the Riemann-Liouville definition (2) inside the integral to 
act on the function f(t): 

x d n f(t) 

°Dzm ee jr^m = j {x _^ a _ n+i dt . (5) 

c 

So, for a < 0, n = 1 c = 0, we have the following expression: 

c n a k _ r ( fc ) kr k-a 

which coincides with RLD D"x k . In general, RLD does not coincide with CD [1]: 
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RLD and CD are not defined for power function x k with k = — 1 due to divergency of the integrals 
in the upper limit. The derivative of 1/x can be obtained as the Weyl derivative defined as 

oo 



T(n -a) J dt 
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which provides D a x~ 1 = — a:~( 1+0! 'r(l + a) for a < l(n = 1). 

The left and the right CD are defined similarly to those for RLD, i.e., due to replacing the 
limits c — > x x — > d in definition (5) consequently. 

One needs to be aware that according to the formulas of addition of orders, the following 
holds([l], p.161): 

n , \—a —j 

D^Dif(x) = D^f(x) -Y J D^I{c+)^0—- ) , 

that is D"D@f(x) ^ D x +l3 f(x), if only not all derivatives D^~ J f(c+) at the beginning of the 
interval are equal to zero. That is why D"D"f{x) ^ D 2a f(x) in the general case. Generalizing 
the Laplace operator in the equation for Newtonian gravitational potential, the author of [5] 
wrongly doubles the order of the repeated fractional derivative. The authors of [11] have avoided 
this mistake, having written down the Laplacian A" as: 

A«u -i-D-fr^tO + r2(a + 1) ^(sin" A + l!^±ll^ 
r 2a Ur {r Ur U) + r 2 « sin" 8 88 ^ 88 ' + r 2 » sin 2 " 8 d<t> ' 

However their research did not concern cosmological problems. 

Let us note one more property of the fractional derivative expressed in modification of the 
Leibniz rule ([l],p.l62): 

oo 

D a x [f(x)g(x)} = ]T m 2_ k + l) D a x D a x - k f{x)D k x g{x), (6) 

which becomes the usual rule as a = n. It can be represented as the integral over the order of 
fractional derivative: 

oo 

d*. [/(«)*(*)] = / r(^ + i)r(al\-M) JrM/(a)g ' g(g) ^ 

— oo 

Later these rules of fractional differentiation will give us an essential modification to the cosmo- 
logical models with fractional derivatives. 

At last we want to note, that RLD can be expressed with the help of the integer derivative at 
the initial point of the interval and CD ([1], p. 163) due to definition (5): 

n / \k — a 

D *fW - E r(i + fe - tt ) /(fc)(gf) + CD * f{x) ■ 



3 Cosmological models with fractional derivatives 

Presumably, for the first time such models within LSM were considered in [2]. To avoid a 
conflict between the occurrence of fractional derivatives in the Friedmann equations and classical 
definition of a tensor in the Einstein equation, based on the integer order derivative in tensor law 
of transformation, the quoted author has repeated the well known derivation of the Friedmann 
equations for a dust from the classical approach (see [3], [12]) and than has replaced all integer 
derivatives with its fractional analog. For the spatially flat Universe, the Friedmann equations 
with a cosmological term are written down in [2] as follows: 

(D?a(t)) 2 = {A x (Gp) a +B{Kc 2 ) a )a 2 , 

(7) 

D?{D?a{t)) = -{A 2 (G P T - B{Kc 2 r)a, 
where a(t) is a scale factor of the Friedmann-Robertson- Walker (FRW) line element, 

ds 2 = dt 2 - a 2 (t)(dr 2 + £ 2 (r)dn 2 ), (8) 
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where £(r) = sin r, r, sinh r for the sign of space curvature k = +1,0,-1, consequently. The 
occurrence of a - degrees in the right-hand side of equations (7) is caused, probably, by dimensional 
reasons. Then the author of the quoted work obtained the following solution to this set in a static 
case a = ao = constant: 

The conclusion is made, that the density of matter and cosmological constant decrease as 1/t 2 if 
G and speed of light in vacuum c remains constant, but if the density of matter and A remain 
constant, then the following holds: G <~ 1/t 2 and c ~ 1/t. The solution of equations (7), mentioned 
as an illustration of the method of [2], in the form a(t) = aoEi ia (Ct a ), where 



t k 



is the so-called Mittag-Leffler two-parametric function [1], is not confirmed by any calculations. 

It is interesting to note that for k = and in the absence of matter and the cosmological 
constant, that is if p = and A = 0, equations (7) are reduced to 

D?a(t)=0, D?(D?a(t))=0, (9) 

which for a = 1 gives the obvious result: a(t) = a = constant, and interval (8) is reduced to those 
of the Minkowski space. If a ^ 1, then the substitution of the zero constant (that is the derivative 
Dfa{t)) from the first equation of (9) into the second one results in identity irrespectively to the 
definition of fractional derivative: RLD or CD. The solution of the first equation in (9) for CD, as 
well as in the case of integer derivative, is equal to constant, but for RLD its partial solution is 
zero, and the general solution depends on time: a(t) <~ i" -1 for < a < 1 (see [1], p. 216). This 
circumstance and also the fact, that initial conditions for equations with CD should be expressed 
by means of integer derivative, instead of fractional one, as in the case of RLD, frequently decline 
a choice of definition for the benefit of CD. In general considerations, we shall not specify a type 
of fractional derivative so long as it will be possible. 

Let us notice that in [5] the Universe dynamics equations are written down in other form, 
namely: 

3{k+(D?a) 2 ] = K pa 2 , a s D?p = D?[a 3 (p + p)}, (10) 

where we use p for the energy density and p for the pressure. The second equation represents 
the energy-momentum conservation law (T^ — 0) for the perfect fluid in space-time (8) with the 
integer-order derivatives replaced by its fractional analogues. Assuming further k — and the 
power-law dependence of a,p and p on the time: a = Ct n , p = At m , p = Bt r , with the help of 
formulas (3) and (10) we can find: 

. R 3 L(n + 1) 2 

m = r = -2a, B - 



k T(n + 1 — a) 2 ' 

and then from the equation of state p = (7 — l)p and the second equation in (10), the following 
equation for degree n in the law of evolution of the scale factor is obtained: 

(7 - 1) L(l - 2a) _ L(3n - 2a + 1) 
j r(l -3a) ~ r(3ra - 3a + 1) ' 

In the same spirit of naive approach, as it was named by the author of [5], the article [6] is 
prepared. In it, the Riemann curvature tensor and, as a consequence, the Einstein tensor are 
defined by the unusual Christoffel symbols containing fractional (of order < a < 1) derivatives 
of metrics coefficients, 

Kx(*) - + dl9 P „ + 9>a), (11) 

where d" is a fractional derivative (2) or (5) with respect to x v . The quoted author writes down 
the Einstein equation, 

. . 1 . . 8irG . . 
R^{a) ~ -g^ v R(a.) = —^-T^a), 
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and the geodetic equation, 

d?x» dx»dx x 

but docs not give any solutions, and only demonstrates that linearized equations in the limit 
a — > 1 are reduced to the usual equations with the integer derivatives for the gravitational waves 
and the Newtonian potential, that was expected from the very beginning due to the definition of a 
fractional derivative. Let us note, that an attempt to prove such fractional derivative replacement 
in Christoffel symbols, such as (11), is undertaken in [5] on the basis of fractional-differential 
geometry, that would be related to FSM formalism. However, the quoted author also mentioned 
that it is so far not clear what such geometry should be. However, in [5] an attempt of construction 
of fractional geometry with the help of fractional coordinate transformations dx % — D^x l dx^ is 
undertaken for the fiat two-dimensional space. 

There were much more advanced and proved results concerning FSM formalism and stated 
by S. Vacaru [13] (see the bibliography therein). In those works, the results of construction of 
the fractional theory of gravitation for the space - time of fractional (not integer) dimension are 
obtained. The author sees one of the simplest motivation for application of fractional differential 
calculus in the theory of gravitation in an opportunity to avoid singularities of the curvature 
tensor of physical meaning due to the completely different geometrical and physical solutions of 
the fundamental equations. Besides, it is noted that models of fractional order are more adapted to 
the description of processes with memory, branching and hereditarity, than those of integer order. 
The result of the application of the method developed by the author of nonholonomic deformations 
to cosmology was the construction of new classes of cosmological models [14]. However, the last 
does not concern to the models under consideration, and it will not be considered here. 

Let us now mention papers [3, 4], where the approach to the dynamical field theories in general 
and to the theory of gravitation is developed on the basis of the variational principle, formulated 
by the author, for the action of fractional order (FALVA). In this approach concerning ISA, the 
integral of action 5l[?] for the Lagrangian L(t, q(r), q(r)) is written down as the fractional integral 
(1): 

t 

SLfa] = f^r / L^qiiT^qiiTW-ry^dT, (12) 

r(a) J 

to 

being at fixed t the Sticltjes integral with integrating function gt{r) = -^n ^_ — y[* a — — T )"L 
having the following scale property: 

9nt{/J-r) = fJ, a g t (T), fi > 0. 
Then (r) satisfies the fractional (or modified) Euler-Lagrange equation: 



dL d fdL\ 1 - a dL 
dqi dr V d<ji J t — T d<ji 



= F , i = l,2,...,n; re(0,i), 



where a dot above the appropriate function stands for the first derivative with respect to time r, 
F 1 is the modified decaying force of " friction" , that is the general expression for non-conservative 
force. In the article [15], time r is treated as the intrinsic (proper) time, and t is the observer time. 

The author of [3, 4], [16] states that at r — > oo we have = 0, and provides some examples of 
application FALVA to the Riedmann geometry and perturbed cosmological models. Every time 
the quoted author does not apply FALVA directly to the gravitational action S , c[a(t)], where a(t) 
is the scale factor in the FRW model (8), but tries to take into account influence of the fractional 
order action (12) on the Friedmann equations through perturbed (and time-dependent) classical 
gravitational constant. Considering Lagrangian L — gij(x, x)x l x : ' , the author obtains the modified 
geodesic equation: 

y ■ '■^■r- I \ 0. (13) 

where T l - k are the standard Christoffel symbols, and T = t — r. The second term here is interpreted 
as a dissipative force, which infinitely increases as r — > t for a ^ 1 and under condition of fixing 
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future time t. Considering equation (13) in the Newtonian approximation, the author established 
the time variation of Newton's gravitational constant (i.e., tried to realize the Dirac hypothesis [17]) 

3(1 Ct) CL 

and introduced the perturbation of the gravitational constant AG = — - — — . Only then, the 

Airpl a 

effective gravitational constant G e ff = G+AG is substituted to the standard Fricdmann equations: 



1,0 , x SlvG e ff A 

+fc ) = -r^+3 

a 47rG e / / , _ , A 



-(a 2 + k) = — ^-p+-, (14) 



(P + 3 P ) + -, (15) 

u o o 

where the cosmological constant A equals zero or A = (f3/t)(a/a), as it is made in [18], where 
G e ff = G = const. Then in [18], [19] the solutions of equations (14), (15) are investigated 
for various equations of state (p — jp), which could make sense if these equations or the way one 
obtains them were enough justified. It would be natural to expect of the author of FALVA the direct 
applications of FALVA to construction of the modified Fricdmann equations from the fractional 
functional Sc[a(t)] that leads to other results, as it will be shown below, without necessity to 
use equation (13) for the mentioned above interpretation of the extra terms in modified equation. 
Moreover, with the help of canonical parameter s = gt{r), equation (13) is reduced to: 

x l +r jk x j x k = 0, (16) 

where the over dots stand for derivatives with respects to s. Actually, the last means that equations 
(13) could be obtained without applying FALVA but simply by replacement of the parameter 
s = g t {r) in equation (16). By the way, it is possible to address precisely the same remark to [20], 
where the modification of cosmology is undertaken on the basis of the periodic weight function 
<7t(r) in action: 

t 

S L [ qi } = J L(T, qi (T),qi(T))exp(-xsm(pT))dT, (17) 

to 

which is named fractional in [20], though such name could be applied to functional (12) only, 
meaning its origination from fractional integral (1). As it follows from (17), in this case the weight 

function is defined by ^ - = exp(— \ sin(/3r)), and replacement s — gt(r) in (16) is immediately 
resulted in equation obtained in [20] by variation of (17): 

x l - Px cos(/3r)x 4 + T) k x>x k = 0, (18) 

again with treatment of the second term in (18) as the perturbation of gravitational constant AG = 
3x(3 cost (3t)H a . . 

, where the Hubble parameter H = — . Certainly, it would be possible to experiment 

47T/3 a 

further with various weight functions to solve those or other problems of the gravitational theory 
and cosmology but it is reasonable before all to return to the framework stated at the beginning 
of our paper and to apply FALVA directly to the gravitational field of the Universe. By the way, 
it was possible to understand the purpose that El-Nabulsi declared while formulating FALVA (see, 
for example, [21], [22] just so). 



4 Fractional derivative cosmology of scalar field 

First we consider the naive (or LSM) approach to the fractional derivative cosmological models 
of a scalar field. For what follows, it is useful to reproduce the derivation of the gravitation and 
scalar field equations from the variational principle for the Einstein-Hilbert action inspired by 
ADM formalism in cosmology (see, e.g. , [23]). The Einstein-Hilbert action-like functional for 
FRW model of the Universe 

ds 2 = N{tfdt 2 - a 2 (t)(dr 2 + Z 2 (r)dn 2 ), (19) 
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where N is a laps function, filled with a real homogeneous scalar field <f>(t), is as follows: 

S EH = S G + S Sc = ^j dtN {-^ +ka-±pj+J dtNa 3 - , (20) 

where V((f>) is a potential of the field. By variation over a(t), (j>(t) and N(t) (with the subsequent 
choice of the gauge N = 1) in the action (20), one obtains the following standard Friedmann and 
scalar field equations: 



'a a 2 k 
2- + — + — 

a a z a z 

a?_ k_ 

a 2 a 2 

On the other hand, equation (23) could be derived with preliminary gauge N = 1 in (f 9), that is 
with proper time and representation of the FRW space-time interval as (8), proceeding from time 
invariancy of the action (20), as it was made in [10]. As a matter of fact, both approaches are 

dL 

equivalent, because (23) was derived from the Euler-Lagrange equation for N(t), that is = 0, 

which just expresses noted time scale invariancy. Besides, instead of equation (22) one frequently 
uses the following equation: 

- = -— (f-V(<f>)) + -, (24) 

which turns out from equations (22), (23). 

Considering LSM, it is necessary to make substitution of fractional derivative of the scale factor 
and the scalar field instead of integer derivatives in equations (21), (23) and (24). Therefore the 



required set of equations becomes as follows: 

^(^+3(^)^+^=0, (25) 

{D> f + k = ^ Q m f + a 2 + A a 2 ; (26) 

Df(Dfa) = ((D^f V(4>)) a+^a. (27) 



While among three equations (21) - (23) of standard cosmology only two equations are inde- 
pendent, and the third one can be derived as the differential consequence of two others, due to 
the modified Leibniz rule for the fractional derivative (6), all equations of (25) - (27) are generally 
independent. The last means that one has either to solve equations (25) - (27) for a(t),(j)(t) and 
V(4>(t)) with constant G and A or to admit the dependence of G and/or A on time when the 
potential V{<j>) is given. This circumstance allows us to doubt the acceptability of models offered 
earlier in the frameworks of the naive approach, which are submitted by the modified equations 
(7) and (10). 

The solution of nonlinear fractional equations (25) - (27) is rather problematic. Even their 
classical prototype (21) - (23) does not always have exact solutions. Unfortunately, the theory of 
equations in fractional derivatives is incomparably more difficult and less advanced in comparison 
with the theory of integer-order differential equations (see, e.g., [1], [24]). Nevertheless, we would 
like to give an example of an exact solution to (25) - (27) demonstrating the fact of existence 
of some solutions. Let us consider the spatially flat model of the Universe without cosmological 
constant (k = A = 0). We assume the dependence of the scale factor, scalar field, and potential on 
time as follows : a = a t", <p — <f>ot m , V(<f>{t)) — Vot r . Due to the modified rules of differentiation 
(3) for the given functions, it follows from (25) - (27) that: 



> • dV(d>) n 
' + 3-0+— ^ =0, 

a acp 



-SnG^j-V^j +A, 



(21) 
(22) 

(23) 
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VW = V *2f*>-(±\ " • (31) 



m = n, r = 2(n — a) , (28) 
n — 2a 

4> n 

The power n in the dependence of scale factor on time is connected to the order of fractional 
derivative a by the relationship: 

in which both parameters appear in the arguments of To/-function, that essentially complicates 
further general analysis of the model. Nevertheless, as the scalar field in equation (29) is real, the 

/ 6 \ 

relationship (32) is possible if n € I —a, a I , where < a < 1. 

Now we consider the modification of the main equations of the classical model (21) - (23), 
which this set undergoes at the intermediate approach (or ISA). For this purpose, we replace the 
derivatives over time in the action (20) with its fractional (of order a) analogous: 



TEH 



j L EH dt = J 



dtN 



(33) 



8nG V N 2 3 J \ 2N 2 

Variational problem with fractional derivatives for functions qj(t) in the action S[qj](t) = 
L(qj (t), c D^qj (t), tD^qj (t))dt on the interval [c, d] yields the modified Euler-Lagrange equations 
[1],[25]: 

dL . „„ / OL \ « ( 8L 



a 

I 



% + ' c nw^)j + ^(^J =a ,34) 

In our case qj(t) = (j)(t),N(t) and a(t). Therefore for Leh from (33) and (34), we arrive at the 
following basic equations of the model: 

t D> 3 D^)-a 3 ^M=0, (35) 

(D?af + k = ^ Q (D?<f>) 2 + V(d>fj a 2 + |a 2 , (36) 

2 t D a (a Dfa) + (Dfa) 2 -k = 8nG ((D?<t>) 2 - V(<f>fj a 2 - Aa 2 , (37) 

where we denote Df = Df and = t D a . If in the action S[qj](t) the derivatives are CDs, then 
(34) remains practically without changes, with the only exception of the derivatives of Lagrangian 
in the brackets are CDs [26]. Hence if the derivatives Df a and D?<j> in (33) are determined as CDs 
they also should be understood as CDs in equations (35), (37) but the repeated right derivatives 
(see definition (4)) in the first terms of (35) and (37) do not undergo redefinition. Let us note 

d_ 

di' 

tD 1 = — — , the set of the fractional equations (35 -37) coincides with the classical one (21 -23). 



also, that in the case of integer order derivative (a = 1) due to the obvious expressions D\ = 



dt 1 
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5 Cosmological models of scalar field with fractional action 



We consider now the cosmological model of a scalar field, which follows from the variational 

principle for the fractional action (12). This approach is in the framework of the intermediate 

approach (ISA) and frequently, as it was mentioned above, is referred to FALVA. Here, it is 

necessary to pay attention to the following peculiar properties of computation to avoid mistakes 

in obtaining the master equations of the model under consideration. The first term in the action 

(20) is obtained as the result of integration by parts of the first term in Ry'—g, which depends 

on the second derivative d. Actually this integration removes the derivative of the laps function 

N from the action. It is a well known procedure. But now, if one makes simply fractional, 

as (12), generalization of the Einstein-Hilbcrt action on the basis of expression (20), the result 

will be incorrect, and some terms will be omitted. Due to the stated above reason, such terms are 

absent in all quoted articles by El-Nabulsi, unfortunately also containing several other inaccuracies. 

t 

Therefore, we use the modified Einstcin-Hilbert action S% H = J L EH dr as the following fractional 

o 

integral: 



t 



T(a) 



o 



a 2 d ad 2 a 2 dN , Aa 3 \ , / ecp 2 



8irG \ N 2 N 2 A 3 3 / \ 2N 2 



(t-T^dT , 



(38) 

where all functions in L^ H depend on the intrinsic time r, and e = +1,-1 for the usual and 
phantom scalar fields respectively. Varying the action (38) over qi — cf),a and N with the subsequent 
choice of the gauge N — 1, we obtain the following Euler-Poisson equations, 



8L% H d (dL% H \ , d 2 (dL EH 



dqi dr \ dqi ) dr 2 



= 



for our model [27]: 



•• /d 1- a\ ■ dVU) „ 

a 1 — a f d\ (1 — a)(2 — a) 8ttG ( Tr , ,.\ A ., , 

a + —{a) + 2t 2 = ^ - M + 3 ' (40) 



2 



I- a f d\ k 8irG / S 2 „, ,A A 



a) + ^UJ + ^ = — v e t + y( ^ + 3' (41) 

where and below the dots above functions designate the derivatives of the appropriate order with 
respect to time t obtained by t — r = T — > t. Precisely the same equations one can obtain from 
generalization of the Euler-Lagrange equations in FALVA for the Lagrangian with higher (second) 
derivative (see [28]) or from the previous integration by parts as mentioned above. In the last case, 
one has to take into account the weight function (t — r) Q_1 in the action integral (38) requiring 
the limit lim t _ s . (o 2 a/i) to be finite. 

It is easy to verify that the set of equations (39) - (41) is represented by three independent 
equations, against two independent equations in the classical case (21) - (23), that is connected with 
violation of the energy - momentum conservation law T^ v v = 0. Therefore in the set of equations 
(39) - (41) for three unknown functions (a(t), <j)(t) and V(t)), potential is not an arbitrary function 
in general. One can rewrite equations (39) - (41) in terms of effective energy density p(t) and 
pressure p(t), taking into account the well known expressions: 

p = + p = - ■ ( 42 ) 

As a result, it follows from (39)-(41) and (42) that: 

P + 3^ + I^W,) = 0, (43) 
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a l-afa\ (1 - a)(2 - a) 4ttG , „. A ,,,, 

a + ^rU) + - i = -— ^+^)+3» (44) 



2 

d\ 1 — a / a\ k 8irG A 



,a J t \a J a 2 3 3 

It is easy to integrate equation (43) for the perfect fluid with equation of state p = -fp: 

= Po 

P 3(l + 7 )i(l+7)(l-«) ' 



(46) 



The interesting fact is that the cosmological A-term in the action (38) can depend on time, but in 
equations (39) - (41), and also in (43) - (45), the only change appears in the dependence A = A(t). 
Therefore the expression (46) is also valid. 

Let us provide an example of an exact solution for the flat model (k — 0) and for the quasi- 
vacuum state of matter: 7 = — 1. From (46) it follows that p{t) = po = constant, as well as in 
the standard model. Then, the remaining equations of (43) - (45) for the Hubble parameter and 
A-term can be copied as follows: 



H-l^H+ ^- a ^- a) = 0, (47) 
2t 2t 2 y ' 

H 2 + — - — H = —po + - . (48) 

From equation (47), it is easy to find that the Hubble parameter varies with time as 

1 — a 

H=^ + Hot 2 , (49) 

where C a = — — a ) ^ jj q j g a cons t an t f integration. From (49), the following dependence 
(3 -a) 

of the scale factor on time t appears: 

/ 3 - a' 

a = aot Ca cxp ^^Hot~^2~ \ , (50) 

while the cosmological A-term changes with time as 

a o H 2A- a , o n (l-<*)(7-3tt) .~ i 3(1 - a) 2 (2 - a) (5 - 2 a ) 1 

A = 3H t +3H Q ^—j t 2 + ^—-- 2 ^-SnGpo. (51) 

It is obvious that at the proceeding to the standard model in the limit a — > 1, the obtained 
solutions (49) - (51) reduce to the known exponential expansion of the Universe : a = a e^ ^, H = 
H , A = 3ff 2 - 87rGp - 

We consider now the dynamics of the flat model of the Universe (k = 0) filled by a scalar field. 
It is convenient to rewrite equations (39) - (41) in terms of the Hubble parameter in the following 
form: 

A-i^H + (1 -^ 2 -°> =-W. (53) 

It is easy to prove that the given set of the independent equations can contain some arbitrariness in 
a choice of unknown functions, for instance H(t) or V(t), only if the cosmological A-term depends 
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on time. However, it is possible to proceed from some dependence A(t). Let us consider one 
example of an exact solution, assuming that the Hubble parameter and the scalar field depend on 
time as follows: 

H=j, <fi=j, (55) 

where C, B are the constants to be defined. We take into account definitions (42) and the equation 
of state p = 7/9 for convenience of representation of solutions in terms of physically clear parameters 
and fractional order of the action a. Using substitution (55) in (52) - (54) and considering (42), 
we obtain the following result: 



ao t 



1-7 
1 + 7 



+ 



a 



where the following notation is used: 



A = A(a, 7 )=± 



■]ni + , 



24ttG 



a 



1+7 



V{4>) 



2e 



e 

2A2 



3 — a 



1-7 
1 + 7 



2a 



(56) 



-1/2 



C = I" 



1+7 



It is easy to see that the baratropic index 7 is not arbitrary but it should satisfy the following 
inequality: 



1-7 
1 + 7 



> 2e 



3 — a 



2a 



and also the inequality which follows from a weak power condition p > 0: 7 > — 1. Thus the 
A-term depends on time t as follows: A = A /t 2 , where 



A = (1 - a) 



(1-7) 2 c(l-7) 



6(1 + 7 ) 2 6(1 + 7) 



(9 - 4a) - 1 



The above mentioned solutions demonstrate the existence of exact solutions to the model, which 
can be of interest in aspect of modelling of processes occurring in the Universe. Let us investigate 
more generally the properties of the masters equations for the spatially flat cosmological model 
(52) - (54). After simple manipulation, this set of equations can be reorganized as follows: 



H + 3H 2 



AnGecf, 2 = 3H 2 



2(4 - a) „ (l-a)(2-a) 



-H 



t 2 



3(5 - a) „ 3(l-a)(2-a) 
-H 



2t 



t_A 

l-a ' 
tA 



2t 2 



1 



2t 



2t 2 



(57) 
(58) 
(59) 



The term in the right-hand side of equation (57) is obtained by dividing the both sides of 

1 — a 

this equality by (1 — a). Therefore equation (57) is valid for all a G (0; 1) but in the classical limit 
a — > 1 it simply means A = 0, i.e. A = constant. The latter is the only consequence of a = 1, and 
among three equations (57) - (59) it is necessary to solve only two of them with A = 0, that allows 
us to define one of the function among Hit), <j>{t) or V((f>) arbitrarily. 

The situation is essentially different for a ^ 1. Now we have three independent equations for 
four functions: H(t), <t>(t), V((j>) and A(t), and specification of A(t) allows us to find the dependence 
Hit) from equation (57). After that, it is possible to find cp{t) and U(£) from equations (58) and 
(59). It is important that the evolutionary, i.e. containing derivative over time, equation (57) for 
the Hubble parameter experiences the influence of the parameter A(t) only. Thus, the behaviour 
of the field and its potential becomes secondary, not determining an expansion dynamics of the 
Universe directly, and can be simply found from the equations (58), (59). In some sense, the 
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scalar field here plays a role of a latent parameter. Indeed, the system behaviour is determined 
by the field behaviour but the Hubble parameter and the scale factor are determined not directly 
by the scalar field but only through the A(i)-term. Using some phcnomcnological expression for 
the cosmological parameter, it is possible to find the dependence of expansion on time, and then 
the scalar field that provided it. There is a certain sense in that, as we know practically nothing 
about a scalar field that evoked cosmological inflation and the present accelerated expansion of the 
Universe. 

We consider one example of the mentioned approach, assuming A = constant. In this case, 
equation (57) contains the only free parameter: a, and can be easily solved, that gives the Hubble 
parameter as: 

1 / 1 - Cn f w a \ 

H= — \9-2a-w a —tjt- , w a = \/16a 2 - 72a + 105 , c > 0. (60) 

6i \ 1 + cot a J 

One can notice, that w a > for a € (0; 1). Integrating equation (60), we obtain the following 
expression for the scale factor: 

1/3 



a = a {t( 9 - 2a - ™«V 2 (1 + c Q t w ») } . (61) 



Substitution of the Hubble parameter from equation (60) in (58), (59) gives us expressions for <j)(t) 
and V(t). Due to the definition (42) and equations ( 58), (59) the energy density and the pressure 
of the scalar field (A = 0) are as follows: 

1 ( q 9 \-c^ \ ( l-cot w ° 
9 — 2a — w n Tn — 15 — 8a 



96^Gi 2 V" ^ a l + c t w -J\ a l + c t w ° 



P = 



96nGt 2 



1 — c t Wa \ / 1 — c t 1 ~ 

9 -2a- w a 10a - 27 - w a -^r- ) - 36(1 - a) (2 - a) 

l + c t Wa J \ a l + c t w <*J 1 A ' 

From these expressions, it follows that the barotropic index 7 = p/p is not constant, and the 
equation of state changes with time. It is possible to show that at the initial moment of time, 
7(0) strongly depends on the fractional order a, but at t — > 00 it practically loses this dependence, 
slowly approaching — 5/7 as a — >• 1. 

It is interesting that if we assume the cosmological term decays as A = ^H, where /3 is a 
positive constant, then equation (ref 56) is reduced into the following one: 

(1 - a - fS)H + 3(1 - a)H> - [2(4 = a){ \ - « ] - ® H - ^ = a ™ - ^ = , 

differing from the same equation in the case ,8 = 0, considered above, only by the constant multi- 
pliers. Therefore its solution is structurally similar to (60), (61). 



6 Conclusion 

In our article, we have given the critical analysis of the main results and approaches known by 
now that are directed towards the application in cosmology of the ideas of fractional differentiation 
and integration. At the same time, we offer the unified consecutive approach to such a problem on 
the basis of the variational principle for the Einstein-Hilbert action which either includes fractional 
derivatives of the scale factor and the scalar field in the Lagrangian, or becomes fractional one. 
The exact solutions of the modified Fricdmann and scalar field equations essentially differ from 
the ones known earlier that could be useful in aspect of modifications undertaken in modern 
cosmology in view of new observational puzzles [29]. We hope that further investigation of the 
obtained equations and its solutions will allow us at last to answer the question: in what degree 
the ideas of the fractional differential calculus are productive in cosmology. 

The author is grateful to Prof. V.V. Uchaikin for the offer of the topic and the valuable 
discussions. 
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